Radiation from an accelerated quark via AdS/CFT 
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In this paper we investigate radiation by an accelerated quark in a strongly coupled gauge theory 
via AdS/CFT correspondence. According to AdS/CFT dictionary, we can read off energy density or 
energy flux of the radiation from asymptotic gravitational field in AdS bulk sourced by an accelerated 
string trailing behind the quark. In the case of an oscillating quark with frequency fi, we show that 
the time averaged energy density is asymptotically isotropic and it falls off as (<?ymA) 1//2 Q 4 /R 2 
, with distance R from the source. In a toy model of a scattered quark by an external field, we 

simply estimate Poynting vector by the bremsstrahlung radiation and show that the energy flux is 
anisotropic outgoing radiation. Based on these investigations, we discuss the properties of strongly 
coupled gauge theory radiation in comparison with electromagnetic radiation. 

<!■ 

PACS numbers: 11.25.Tq, 04.20.-q 

" I I. INTRODUCTION 

It is quite interesting to investigate the mechanism of energy dissipation by an external quark 
moving through gauge theory plasma from both theoretical and experimental point of view. Recently, 
a strongly coupled gauge theory plasma is believed to be formed at RHIC [H, ■ A heavy quark 
£N1 ■ moving through the plasma is considered to lose the energy via two-body collisions with thermal 

quarks or gluon emission. Although it is detectable through jet-quenching effect at RHIC, this 
energy loss rate cannot be calculated analytically beyond the perturbative QCD. 

Gauge/gravity duality [!, 0, Hj] shed some light on the investigation of the strongly coupled gauge 
theories. This duality states that the expectation value of the energy-momentum tensor ( ) of 
thermal Af — 4 super- Yang-Mills (SYM) theory at large 't Hooft coupling gy M -/V and a large number 
of colors N can be read off from the perturbed metric of anti-de Sitter (AdS) black hole. For example, 
we can obtain transport coefficients such as shear viscosity by calculating quasi-normal mode of AdS 
black hole. Many calculations showed that the ratio of viscosity divided by entropy density is equal 
to 1/47T for a large class of gauge theories described by gravity duals (See, for example, Ref. 0,0] 
for all references). This universal value is in good agreement with the RHIC data, which suggests 
that it is a universal feature for all strongly coupled gauge theories, including QCD. This agreement 
' leads us to seek other universal properties of the strongly coupled plasmas. Even though there is 

no universality, we would obtain some insights which are not readily obtained from perturbative 
techniques and lattice simulations. 

In gauge/gravity duality, a heavy quark corresponds to an open string which stretches from the 
boundary of AdS spacetime to the horizon of the black hole. Under the test string approximation, 
a large number of studies have been made on estimating a jet-quenching parameter or energy loss 
rate of the quark moving through the plasma via gauge/gravity duality [a. M. ITol. [TP. [TH. [l3l. UA. [l5| . 
Beyond the test string approximation, we need to take into account the back-reaction by the radiation 
from the string onto the original AdS geometry, as pointed out in (l5j . So, a good place to start is to 
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investigate the properties of radiation emitted by the string. In a series of papers [ltil. H?! [HI [l9l. [20| . 
the energy density and energy flux associated with the quark with constant velocity were calculated 
from the metric perturbation by the string trailing behind the quark. In the case of supersonic 
motion, the Mach cone structure, which is characteristic of Cherenkov radiation, was also displayed. 

In this paper we shall investigate the properties of radiation from an accelerated quark via 
gauge/gravity duality. Since the quarks in the plasma would be scattered by one another, the 
bremsstrahlung radiation by the scattered quark would bring about another mechanism of energy 
dissipation of the quark moving through the plasma. Callan and Giiijosa [2l| first investigated time- 
dependent fields produced by an oscillating quark via gauge/gravity duality and showed that the one 
point function ( TrF 2 ) in the boundary theory decays faster than 1/R 2 . This does not immediately 
mean that there is no radiation emitted by the quark. As mentioned in [2ll [22j , we can explain the 
reason for the curious behavior of the function by analogy with electromagnctism as follows: F 2 is 
expressed by the electric field E and magnetic field B as F 2 = E 2 — B 2 . Since E 2 = B 2 ~ 1/R 2 
for radiation field, the leading radiation parts are cancelled in F 2 . So, we show the existence of the 
radiation by dealing with energy-momentum tensor of the gauge theory generated by an accelerated 
quark. 

As a first step, we calculate energy density or energy flux of radiation from an accelerated quark at 
zero tem per ature by solving master equations of linearized Einstein equations sourced by the string 
given in [l9( . In section [IT] we briefly review the master equation and the general method to solve 
the equation developed in . In section [HI] we calculate the time averaged energy density ( T 00 ) 
of radiation in the case of an oscillating quark with frequency Q. It is shown that the energy density 
is asymptotically isotropic and it falls off as y/ gy M N f2 4 /i? 2 . In section [TV] wc simply estimate the 
Poynting vector associated with the bremsstrahlung radiation in a toy model of a scattered quark 
by an external field. It is shown that the energy flux is anisotropic outgoing radiation. Conclusion 
and discussions are devoted in the final section. 



II. MASTER EQUATION FOR A SCALAR MODE 

A quark in M = 4 SYM theory at zero temperature corresponds to a fundamental string through 
the Nambu-Goto action 

S = -^Jd^V=h, (ILl) 

where h is the determinant of the world-sheet metric h a p = g^daX^dfiX 1, '. The background metric 
g is given by the metric of AdS§ spacetime 

ds 2 = a(r) 2 (-dt 2 +dx 2 ) +dr 2 /a(r) 2 , a(r)=r/L, (i = 1,2,3), (II.2) 

where L and a' are related by 't Hooft coupling <?ym-^ as 

LVa' 2 = g 2 M N . (II.3) 

We solve the linearized Einstein equations in the metric sourced by the accelerated string. Since 
we are interested in the energy density or energy flux of the bremsstrahlung radiation, we shall focus 
on the scalar perturbations of the linearized Einstein equations. Following [l9l ]. the scalar modes of 



the string energy-momentum tensor Tu V and the metric perturbations are expressed as 



T ab {t,x,r) = du> 



d 3 k 



T^ c (u,k,r) S(k,x) e 



Tt,i(t,x,r) = a(r) / duo 



d 3 k 



t 6 s (w, k,r) Si(k,x) ] e 



(II.4) 
(II.5) 



(2tt) 3 
f d 3 k 

Tij(t,x,r) = 2a 2 (r) / duo / J77 -^ [ p s (uo,k,r) S(fc, as) + t s (uo, k, r) E tj (k,x) ] e~ wt , (II.6) 



' 'j 
and 



h a b(t,x,r) = du 



(2tt) 3 



d 3 k 

f* b (uo,k 7 r) §(k,x) 



(27T) 



hbi(t,x,r) = o(r) / 



hij(t,x,r) = 2a (r) / cLj 



d 3 fc 
(2^)3 



(27T) 



(II.7) 

3 - [flf(w,fc,r) <% S(fe,cc)+F|(a;,fe,r) S«(fe,x) ] e~ lwt , 



(IL9) 



where a,b = t,r and functions S, Si, and Sjj are defined as 



S(fc,a:) 



Ah-x 



Si(k,x) 



diS(k, x) 



Sij(k,x) = — didj §(fe,x) + - 5jj S(fc,cc) 
Here k is defined by fc := |fe|. 



The diffeomorphism-invariant quantity of f^ b is constructed as [l9j 



/ah - /o6 + 2D (a X b) , 

Defining the master field 1 as 



Xn 



a(r) 



a(r) 



k 2 I ftt _3 IS , a3 ffr 



the master equation is given by 



\ d r+l d r-^[(k 2 -U 2 ) + ^) 



(11.10) 



(11.11) 



(11.12) 



(11.13) 



with a source term of the string 
j(uj,k,r) 



ZL 3 uor 3 



3i 3 



k 2 L 2 uo r 3 T rr (uo, k, r) + 2lo{3L 4 lo 2 r - Ak 2 L 4 r + 48r 3 )r 5 (w, k, r) 
- ikL{ - 3iuo r 3 r r s (w, k, r) + (fc 2 L 4 - 9r 2 )rf (lo, fc, r)} 

2k 2 L 2 Tg{ij 1 k,r) + kLr{kLT J J(u), k,r) - 3t{ s (cj, k,r)} 
+ 6ir 2 oo {9t' s {oj, k, r) + r t" s (lo 7 k, r)} 



(11.14) 



1 We set k\ = 1. 
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where a prime means the derivative with respect to r. 

Eq. pi,13[) can be solved by a Green function Gr(uj, k, r; ro) which obeys 

CG R (u},k,r;r ) =S(r-r ) . (11.15) 

This is constructed by two solutions of the homogeneous equation £$>h = = £$b, where 
satisfies appropriate boundary conditions at the degenerate horizon r = 0, while &b has appropriate 
falloff at the boundary of AdSs . Then the Green function Gr satisfying the boundary conditions at 
both the horizon and infinity is represented by <f>s and <f># as 



G R {u},k,r;r ) 



^ H {r )^B(r)e(r -r ) + ^>H(r)<P B {ro)0(ro -r) , (11.16) 



W[^ H ^ B ](r Q ) 
where W[&h ,&B](ro) is the Wronskian defined by 

W[$ H , $ B ](r) = $ H (r)d r $ B (r) - <S> B (r)d r <5> H (r) . (11.17) 
So, the solution of Eq. (|II.13[) is given by 

$(w,fc,r) = A(r)$ H (r) +B(r)$ B (r) (11.18) 

with 

-^-f *°t^ i(r ° ) ' B < r >— «•»>■ (,u9 » 

As discussed in [2(|, the expectation value (T^„(i, aj) ) of the energy-momentum tensor of jV = 4 
SYM theory can be read off from the asymptotic behavior of 3>(ui, k, r) near AdS$ boundary, 

*(*,*, r) = + + 0(r" 3 ) . (11.20) 

The term P/r in Eq. pi.20p produces singularity of (T^ v (t^x)) supported at the location of the 
quark, while the Q/r 2 term gives the energy density of the radiation as 

(T 00 (oj,k)) = Q(uj,k) . (11.21) 



III. RADIATION FROM AN OSCILLATING STRING 

In this section, we calculate the energy density ( Too ) of Af = 4 SYM radiation from an oscillating 
quark. This can be obtained by solving Eq. pi,13[) for an oscillating open string. 

Now we consider that the string oscillates only in the x\ direction [X<x = X3 = 0). Under the 
static gauge t — t, a = z = L 2 /r, the equation of motion is obtained from Eq. as 

= . (III.l) 




Here a dot and a prime mean the derivatives with respect to t and z, respectively. Assuming that 
the oscillations are very small (\Xi\ <C 1), the linearized equation is simply given by 

2 

X'l - -X[ - X 1 = . (III.2) 

z 
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Under the ansatz Xi(t, z) = sm{Q£)^{z) , it is easily checked that 
£(z) = e {cos(fiz) + fiz sin(fiz)} , 



z = L 2 /r , 



(III.3) 



is a solution of Eq. (|III.2|) . where e is a small parameter and fl is a positive constant, representing 
the angular frequency of the oscillating string. Calculating the energy-momentum tensor t^ v , 



8S 



-9 % 



2ira' 



h al3 (d a X> 1 ){d X 1 ') Six 1 - X 1 ) 8(x 2 )S(x 3 ) 



(III.4) 



one gets the non-zero components up to the second order 



1 



d 3 k 

47ra'a 5 (r) J (2ir) 3 

-l r d 3 k 



Ah-: 



ih-: 



Aira'a(r) J (2ir) 3 

cos(fit) sin(fii) f(r)£'(r) 
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n 

2ira'a(r) 

n 

27ra'a 5 (r) 
1 



2 - 2ifa sin(fit) f(r) + fi 2 cos 2 (Oi) £ 2 (r) 

+ a 4 (r) sin 2 (fit) £ ,2 (r) - k\ sin 2 (Oi) £ 2 (r)) , 

^2 + 2 j fc x sin(fit) £(r) - fl 2 cos 2 (fit) £ 2 (r) 

- a 4 (r) sin 2 (fii) £ ,2 (r) - k\ sin 2 (fii) £ 2 (r) N 
d 3 fc 



_ife-EB 



(2tt) 3 



cos(fti) £(r) -ifei sin(fit) £(r 
sin(m) ^'(r)^l — ifei sin(Ot) £(r)) 



(27)3 
d 3 /fc 



2na'a{r) 1 'J (2tt) 3 

27rQ/ 1 Q5(r) cos 2 (m) £ 2 (r) - a 4 (r) sin 2 (ta) £'(r) 



d 3 k 
J2^f 



(III.5) 



where a dash denotes the derivative with respect to r. Substituting these into Eqs. (|II.4|) - (|II.6|) . the 
components t s , P s , rf , rf, rf r , and rf r are 



s id kit; 

Ana' k a 2 (r) 

Ana' k a 2 (r) 
1 



S(lo + 2fi) - <5(w - 2ft) 



P 



2Ana' a(r) 



S(uj + fi) - 5{uj - Q) + fci £ %) - ki £ 
> 2 / ft 2 



5{w + 20.) + 8{lj - 20) 



, 2 \ 8{uj + 2Ct) + S(lo - 2Ct) 



3fcf_ 

k 2 



1 P" 



T S , = 



8na' a 5 (r) 



(a 4 (r)£' 2 + k\e + fl 2 £_ 2 - 4) 5(u) + h£ 



8(lo + Q) - 8{lo - ft) 



+ {n 2 e-aHr)e-k 2 e) 



8{uj + 2ft) + 6(lu - 2ft) 



id 
8na' a(r) 



£f (8(LJ + 2fl)-8(ui-2Cl)) 



(111.6) 



G 



We are interested in the time integral of the energy density, 



dt (Too(t,x)) = 



dt 



duj (Too(u,x)) e 



2ir (T 0O (u = 0,x)) 



(111.7) 



Owing to the linearity of Eq. pi.l3[) . only the terms with S(w) in Eq. pil.6[) contribute to the time 
integral of the energy density. So, we drop the terms with S(u> ± 0) or 6(u> ± 2f2) in Eq. pil.6p . 
Substituting Eq. (|III.3[) and Eq. (|III.6[) into the source term (|II.14|) . one obtains 



j(w,k,r) = j (k,r) 5(u) H , 

io(fc,r) = e 2 (j± + ^ + cos(2fiz) 



(mi 



where all the coefficients are given by 



ft 
ft 
Si 
iVi 
Aa 



ft 



167ra' 

l 8 n 2 

487ra' 
L 2 Q 3 
87ra' 

n 2 

167ra' 



C 2 = 



L 4 



3fcf_ 
fc 2 



rifcf + 8ft(^ + 3ft) l- 



487ra' 



AA 



(III.9) 



18ft -i + 6ft(fcf - IT) + rfjfcf - 3ft) 



3fcf_ 
fc 2 



3fcf_ 

fc 2 
3fcf_ 

k 2 



S 2 = 
iV 2 = 



L 6 n 

24na' 
L 4 k 2 
487ra' 



ft 



k 2 k\ - 3ft(fc 2 - 3fc 2 ) + 6-p-(fc 2 - 3fc 2 ; 
6ft | + * 2 



3ft' 



L 8 n 



48ira 



7 fc 2 (jfcf + 2ft) 



(111.10) 



The two- independent homogeneous solutions and <&b for w = are simply represented by the 
modified Bcsscl function as 



*B(r) 



h(kL 2 /r) 



®h{t) 



K (kL 2 /r) 



and their Wronskian is simply given by 



W = -1/r 5 



(III.ll) 



(111.12) 



We can obtain the coefficient Q in Eq. pl.20|) through a straightforward but tedious calcula- 
tion (the detailed calculation is given in the Appendix A). The coefficient Q is given by 



Q(w,fc) = Qo(fc) <S(w) + " 
Te 2 ft 



(III. 13) 



Qo(fe) 



3^e 2 



ft 



2L 2 ft (4 + fc 2 /^ 2 ) 1 / 2 2L 6 ft (4 + fc 2 /ft) 5 / 2 L 6 ft (4 + P/ft) 3 / 2 



Si In 



2ft 
~k 



4ft 



+ 



kL N\ 



f 1 
N 3 



s 2 n 



ft (fc 2 +4ft) 3 / 2 



No 



kL 2 (kL 2 ) 3 



(III. 14) 
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The expectation value (Tqq(cj,x) ) in position space x is obtained by the Fourier transformation, 



r°° rfh 

(Too(o;,x)) = J ^ Q(u,k) e ikx . (111.15) 



'-co (2^) 3 

So, in the far asymptotic region \x\ 3> 1, the integration is dominant for k — > 0. In this limit, 

^e 2 / , 3C 3 k 2 C 3 

2 9 

" — Sllnk+ 8LHV S2 

The integral can be easily performed by differentiating appropriately modified versions of the stan- 
dard formulas 

d 3 k e lkx 1 f d 3 k e ikx 1 

3 ' 

Ft 

(2tt) 3 |fe| 3 2tt 2 ' J (2tt) 3 |fc| 4 " 8tt ' (HL17) 

where i? := |a;|. So, Fourier transform of every term except last term in Eq. (|IIL16|) decays faster 
than 1/i? 2 2 . The Fourier transform of the last term becomes 



(2tt) 3 |fc| 27r 2 i? 2 ' J (2tt) 3 |fe| 2 4irR ' 

d 3 k e lkx InR f d 3 k e ikx nR 



^"■•"-SS^'M-^^EF'M- < III18 > 

Substituting this into Eq. (|III.7|l . we obtain 

* 96 ff™ (2^(0)) = 96 ff^ • (HI.19) 

This relation implies that the coefficient of S(uj) in Eq. piI.18|) corresponds to the time averaged 
energy density of the radiation. The energy density ( Too ) can be expressed by the acceleration, 
v~en 2 , as y/^iV v 2 /R 2 . 

Let us compare Eq. (|III.18|) with the times averaged energy density of electromagnetic radiation 
sourced by an oscillating electron along x\ axis with the same amplitude e and frequency f2. In 
terms of charge e of the electron, the power radiated per unit solid angle becomes 

e 1 2q4 „2-2 

5W tf9 "T tfe ' (,,I20) 

where is the angle between x\ axis and x. Then the time averaged energy density is obtained by 
Eq. pil.20p divided by the speed of light. 

So, both of strongly coupled gauge field and electromagnetic radiation densities agree on the 
asymptotic falloff and the acceleration dependence, while the angular distribution and the coupling 
dependence in the energy density are quite different. As seen in the case of a static quark [23| . the 



2 In fc term in Eq. 1)111. 16|l also decays faster than 1 / R 2 . 
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energy density in Eq. piI.18[) is proportional to the square of 't Hooft coupling, yj g YM -/V, which is 
characteristic of strong coupling. We will argue on the angular distribution in Sec [Vj 

To investigate the energy loss rate of the accelerated quark, we need to obtain Poynting vector 
( T 01 ) . Unfortunately, it is technically difficult to explicitly calculate ( T Ql ) because of the com- 
plicated source pi.l4|) for the oscillating string solution (|III.3[) . So, in the next section, we simply 
estimate the Poynting vector ( T 0t ) by using a toy model where a straight string which moves at 
constant velocity and then suddenly stops. 



IV. BREMSSTRAHLUNG RADIATION FROM A SCATTERED QUARK 



In this section, wc shall consider brcmsstrahlung radiation from a string which moves at constant 
velocity v along xi-axis, suddenly stops at the origin and then remains at rest. The Fourier transform 
j(t, fc, r) of the source j(ui, k, r) in Eq. 13|) can be constructed by the source term J(t, k, r; v) of 
the string moving at constant velocity v for — oo < t < oo as 



j(t, fc, r) = J(t, fc, r; v)O(-t) + J(t, fc, r; 0)0 (t) , 



(IV.l) 



where the time when the string stops is set to 0. 

So, using Fourier transform j{w, k, r) of j(t, fc, r), we rewrite the master equation in Eq. 13|) as 



= -j(u,k,r) 



J{k,r;v) J(fe,r;0) 



uj± = uj ± it 



o- — k\V uj + 

where the source term of the string moving at constant velocity v is given by 

= ^ + - 2 ) " ^k\} {- A + ^ (2fc 2 - ^k\)) 

1 1 

l(v) 



(IV.2) 



J(k,r;v) 



27ra' y/Y~x? 



(IV.3) 



In order to construct two independent solutions $s which satisfy the homogeneous equation 
= as well as the boundary conditions at the horizon and infinity, we introduce "radial wave 



number" in the bulk, q 



fc 5 



,1/2 



For the multivalued function q in the complex u-plane, we 



take the branch cut such that q is holomorphic in the upper plane, that is, 



y/uJ 2 - k 2 



q = <i \Jk 2 — , 



k 2 



for uj > k 

for k > ui > —k 

for — k > uj 



(IV.4) 



Then and are given by the Bcsscl functions as 



* B (r) = \jo(qL 2 /r) , 



1 



<f H (r) = -4 1) (qL 2 /r) , 



W = A . (IV.5) 



Following to the procedure in the previous section, we obtain the coefficient Q in Eq. pi.20p as 



L 2 



8a'k 2 



{3v 2 k 2 - k 2 {2 + v 2 )} (fc 2 + 3v 2 k 2 - 3uj 2 ) 2k 2 (k 2 - iuj 2 ) 



3q(tj_ — kxv)\l — v 2 
Supposing that \v\ <C 1, we shall expand Q in terms of v as 
Q = (terms up to v) 



3uj+q 



(IV.6) 



, L 2 2k 2 k 2 (k 2 - 3uj 2 ) + |2fc 4 + 9k 2 uj 2 + 3fc 2 (fc? - 2w 2 ) W 

v —r — , l2 3 — — + o(v 3 ) 

a' 24k 2 uj±q 



(IV.7) 
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The expansion would be valid for \u>\ > \ki v\. So, hereafter, we focus on the high frequency behavior 
of the Poynting vector and we simply replace uj- with uj. 
Since the Poynting vector (T 0i (oj, x) ) is defined by as 



T 0i (uj, x) 



d k kiW 



Qe 



ikx 



d 



crfc uj 



(2tt) 3 k 2 ^ ' dx 
the total instantaneous power radiated by the quark is given by 

P(t) = I (T 0l (t 
Jav 



(2tt) 3 k 2 



Qe 



ikx 



(IV.8) 



(IV.9) 



where rii is a unit normal vector to two-dimensional spacclike surface dV enclosing the quark. The 
first order terms of v in Eq. (|IV. 7|) does not contribute to the total instantaneous power P, since 
(T 0l )ni «-> —(T 0l )ni for xi <-> — x\ for the first order terms. Therefore, we only calculate the 
second order terms of v. 

After a little bit tedious calculation (in detail, see Appendix |B|) . one gets the Poynting vector 
( T° l ) as 



[T oi (uj,x)) = (terms up tow) 



iL 2 v 2 d fI (R) Id 



4a' dx 1 
+0{v 3 ) , 



1 d 2 



3 3^^+2^ /2 ^-^ 



3a; 2 d 2 
~~dx- 



MR) 



(IV.10) 



where R = \x\ and Iq, I\, I2, and I4 are functions of R and uj defined by 



i (R) 

h(R) 
h{R) 



AttR 

i 

AttR 
1 



' h[ 1] {\lu\R) for lu > 

H[ 2 \\uj\R) for uj<0 ' 

f -H^\\uj\R) for lu > 
H {2 \\uj\R) for uj < ' 



AttRuj 



MR 



dz x 



for uj > 
for tj < 



h(R) = —5 



fh(f)df df + d(uj)R + C 2 (to) 



(IV. 11) 



with two arbitrary functions C\ and C2 of uj. 

In the asymptotic region R — > 00, the Hankel functions {\uj\R) behaves as 



H$\\u\R) 



H£\\u>\R)~-i 



2 e*'"^ 



7T y^4R 

2~ e j|w|JJ 



^ 2) (kli?) 



2 e 



— ilcJ I i? 



7T 

2" e -i|w|i? 



7T V 7T 

So, in the frequency w-space, the leading order of ( T 0i (cj, x) ) becomes 

,|3/2 



(T^(uj,R))^- 



iL 2 v 2 ujx l ( 2x 



Uira'R 2 V R 2 



- 1 



2irR 



= 



R 3 / 2 



(IV.12) 



(IV. 13) 
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The factor e lujR shows that the energy flux is an outgoing wave, as we expect. Since this factor 
oscillates heavily with respect to u for large R, we might expect that the slow falloff with R in 
Eq. (|IV. 13|) disappears in position space except the neighborhood of future light cone t ~ R. In fact, 
the first term in Eq. (|IV.10[) . for example, becomes in position space 

M6>(t-i?) R6(t-R) 6(t - R) 

(t 2 - i? 2 ) 5 /2 " (*a _ #2)5/2 - t 4( X _ (fl/ t )2)5/2 ' 

where the step function inevitably appears from the causality condition. It is easily seen that in any 
timelike directions (t > R) the term rapidly decays as t~ A . One can also see that the other terms in 
Eq. (|IV.10p decays similarly. So, the slow falloff with R in Eq. (|iV.13j) disappears in the real space. 3 
The strong divergence on the light cone t ~ R may be due to the simple model in which the string 
suddenly stops. 

In terms of 't Hooft coupling gy M N, Fourier transform of the total instantaneous power P(u>) is 
obtained from Eq. pV.13[) as 

PH = / ( T ^, R))n^^J M SZ#/2 e -fl. (IV.14) 
Jdv 9 » 2lT 

It is interesting to compare Eq. (jlV. 14j) with the known result in the case of electromagnetism. 
In the electromagnetic case, we consider the bremsstrahlung radiation sourced by an electron with 
initial velocity v which rapidly stops in the interval At. In the low velocity limit v — > 0, Fourier 
transform of the total instantaneous power radiated is given by 

P(w) = l{ T°>, r) ) m * e^ R . (IV.15) 

Eq. (|IV.14[) does not include the time interval At since the quark suddenly stops with almost zero 
time interval. So, we must compare Eq. pV.14[) with Eq. PV.15P in the range \kiv\ -C <C l/At. 
We find that both the coupling dependence and the frequency dependence are quite different between 
the strongly coupled gauge theories and electromagnetism. 

V. CONCLUSION AND DISCUSSIONS 

We have examined radiation by an accelerated quark in J\f = 4 SYM theory with large 't Hooft 
coupling via gauge/gravity duality. In the oscillating quark with frequency fi, we have shown that 
the time averaged energy density far away from the quark correctly falls off as f2 4 /i? 2 . The f2 
dependence also agrees with the result in electromagnetism. 

Interestingly, the angular distribution of the time averaged energy density is isotropic. Since Af = 4 
SYM theory is in deconfincd phase, gluons can be emitted. Although the energy density is composed 
of not only gluons, but also other fields, it is natural to consider that the dominant ingredient of 
the energy density is gluon. In the weak coupling limit, the gluon radiation should become dipole 
radiation analogously with electromagnetism. So, one may expect that dipole radiation occurs even 
in J\f = 4 SYM theory. In the strongly coupled gauge theory, however, a gluon emitted by the source 
quark induces production of gluons around it by the self-interaction. As a result, we may expect 
that the energy density is asymptotically averaged over all the angles by the gluons. 



/oo 
dui e~ 
-oo 



3 Since the expansion of Eq. l)IV.6|l by v breaks down near the zero mode, we cannot strictly derive the energy flux 
in position space. From the dimensional analysis, however, we expect the power law decay in time. 
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In the toy model of a scattered quark, we derived the total instantaneous power sourced by 
a quark which moves at constant velocity v, suddenly stops at the origin and then remains at 
rest by calculating Poynting vector (T 0i ) via gauge/gravity duality. In this case, we found that 
both the coupling dependence and the frequency dependence are quite different from the case of 
electromagnctism. The Fourier transform of the total instantaneous power in electromagnetism 
P(ui) is independent of the frequency oj, while P(u>) in J\f = 4 SYM increases with uj. 

Even though the static force in strongly coupled gauge theory is well understood in terms of 
Wilson line probe, the radiation problem is still largely open question. In this paper we have 
investigated properties of radiation at zero temperature as a first step toward the understanding of 
the properties at finite temperature. It is still an open question whether or not universality exists in 
the radiation properties, although we have seen it in the ratio of viscosity divided by entropy density. 
If universality exists, we expect that the investigation at finite temperature may be connected to 
the RHIC data through jet quenching phenomenon. Even though there is no universality in the 
properties of radiation, we would obtain some insight into the understanding of the mechanism of 
energy dissipation by an external quark moving through the thermal plasma. 



We would like to thank M. Natsuume for useful discussions. We are grateful to K. Matsui for 
comments. 



APPENDIX A: CALCULATION OF COEFFICIENT Q FOR OSCILLATING STRING 

Since the coefficient Q in Eq. pi.20p is finite part of B in Eq. (|II.19[) . we give calculation of B in 
this appendix. 

Let us define joc, Jos, an d Jon as 



Acknowledgments 




(A.l) 



The coefficient B$ coming from the partial source jos can be calculated as 



B S (r) 




r 



dr + K (kL 2 /r ) sin(2ttL 2 /r ) 



\r r%) 



POO 

SieV7r£ 2 f2 / dz a z~ 1/2 K (kL 2 z )J 1/2 (2L 2 nz ) 



Jl/r 




+S 2 e 2 VirL 2 n / dz Zq /2 K (kL 2 z ) J 1/2 (2L 2 n Zo ) 



(A.2) 



So, in the large r limit, we obtain the coefficient B$ as 




I/ 4 (k 2 + 4ft 2 ) 3/2 



irnS 2 e 2 



(A.3) 
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with the help of Appendix C. 

The coefficient Bn coming from the partial source joiv can be calculated as 



B N (r) = e 2 [ dr ( N, + ^ + ^f) K (kL 2 /r 
Jo \ r r / 

'kL 2 N x N 2 N 3 y 2 

«'</ [ 

JkL 2 /r 

-e 2 N 1 kL 2 



e 2 N 2 
' kL 2 



y 2 kL 2 (kL 2 ) 3 

K (y)r + r 

. kL 2 /r JkL 2 /r 

- 2 N 3 



y 



dy 



kL 2 /r 



dy K (y) 



{kL 



2 V> 



kL 2 /r 



K Q (y) 

Ki(y) \ 
y J 

dy y 2 K (y) 



= -e 2 N l kL 2 
+e 2 (^kL 
where we used the formulas 



K (y) 



kL 2 /r 



[Ki(y)} 



kL 2 /r 



2,^ 

kL 2 



e 2 N 3 



dy K a (y) , 

kL 2 /r V ftij ) JkL 2 /r 



dy y 2 K ( y ) , (A.4) 



K Q (y)-K 2 (y) = --K 1 (y) 

y 



K (y) + K 2 (y) = -2K[(y) , 



(A.5) 



in the third and fourth lines, respectively. In the large r limit (?' — > oo), we obtain the coefficient 
Bn as 



B N = -e 2 kL Ni 



No 



kL 2 (kL 2 ) 3 



0(r lnr) 



(A.6) 



with the help of formulas J °° K (y)dy = J °° y 2 K ( y )dy = ir/2. 

The coefficient Be coming from the partial source joc can be calculated in the large r limit as 



B c (r) 



r 

€ 2 kL , 



dr [ C x + % + ^ ) K (kL 2 /r Q ) cos(2i 2 O/r ) 



'o 'o 

M 2 _l r< s/2\^ 



2 V O Ji 

3 



It 



dz (C 2 z^ + C 3 zl' z )K 1 (kL 2 z )J 1/2 (2L 2 nz ) 



e 2 C 3 In 



/ rOO POO 
- dz z* /2 K (kL 2 z Q )J 1/2 (2L 2 nz ) + e 2 C 1 L 2 dz 
"il/r JL 2 /r 



K (kz ) 

2 cos(2S2z ) 



ire 2 C 2 



ire 2 C 3 (k 2 - 8n 2 ) 



2L 2 Vk 2 + 4fP 2L e (k 2 + 4ft 2 ) 



5/2 



<?C X L 2 



dz cos(2ftz ) , (A.7) 

L 2 /r Z V 



with the help of Appendix C. The third term in the last line includes divergent terms. So, using 
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Eq. (|C.3|) , we carefully calculate as 

P(r) := e 2 dL 2 f°° dz ^° ( ^ o) cos(2^ ) 

i„ , f°° , cos(at) e 2 C x r f°° dt . . .. . 

= e 2 C 1 rcosb dt H -!-/ {(at + b)si(at + 6) + (at - b)si(at - &)} 

f 00 dt { b 2 b 2 1 

~ e 2 Circos6A (a) — e 2 Cir / <atsi(at) — — cos(ai) — — sin(at) > 

io v* 2 + l I 2 2t J 

/ Z" 00 /<\M ft 2 b 2 f a \ 

= e 2 C ir K (a) cosb + a dx — + —K (a) + — / dx K (x) , (A.8) 
V J a x 2 2a J J 

where a = kL 2 /r and 6 = 2flL 2 /r and the function si(x) is defined by 

, , f°° shit 
si{x) ■=- dt . (A.9) 

In the third line, we used the Taylor expansion of si (a;) around x = at for large r. In the asymptotic 
region, r —* oo, we obtain 

P~0(rlnr). (A.10) 

This logarithmic divergence does not appear in Eq. (|II.20|) because the coefficient A in Eq. (|II.19j) 
also includes the same divergence with an opposite sign. So, P does not contribute to Q(cj, k) in 
Eq. (111201) . 
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APPENDIX B: CALCULATION OF { T 0i ) FOR A SCATTERED STRING 

The calculation of ( T 0t ) can be easily computed by differentiating appropriately modified versions 
of the following Fourier transforms of 1/q and 1/kq: 



Io(R) 



h(R) 



1 



(2tt) 3 
1 



d s k 



{2ir) 2 iR J_ 
1 

(2ir) 2 iR 
1 



di- 



ke 



ikR 



<1 



ke tkR sgn(uj) 
dk — — 
j\ \/uj 2 — k 2 



dk 



8(\k\ - \uj\)ke tkR 



(2ir) 2 iR 
1 

(2w) 2 iR 



muJi(\u\R) - — . 

or /_ 



dk 



iVk 2 - uj 2 
6{\k\ - \uj\)e lkR 



Vk 2 - uj 2 



4itR 
1 

(2^)3 
1 



muJi(\u\R) + ir—N (\u\R) 

j H[ 1] (\lo\R) for lu > 
{H[ 2) (\(j\R) for uj<0 ' 

d 3 k " 



qk 



(2n)HR 
1 

(2n) 2 iR 



ikR 



dk 



sgn(w) 



-\u\ V^ 2 - k 2 



dk 



6(\k\ - \uj\)e lkR 



iVk 2- 



Trsga(w)Jo(\w\R) +niN (\uj\R) } 



-H^\\uj\R) for uj > 



^R " li4 2) (Mi?) 



for w < 



(B.l) 



where R = \x\. So, Fourier transformation of l/k 2 q is easily calculated as 



l 



(2tt) 3 
1 

27T 2 i? 

1 

2^R 



Ak-a 



R Jo 
i 



R 



d A k 

DO 

dk 

R 

df 

df rl\(f) 
\R 



k 2 q 

sin(fci?) 



kq 



dk 



cos(kf) 



AttR\l 



dz x 



Similarly, from the fact that 



dR 2 



R 



(27T)* 



d A k 



Ak-x 



k 4 q 



for uj > 
for cj < 



= -RI 2 {R) 



(B.2) 



(B.3) 



15 



Fourier transformation of l/k^q can be represented by two integration constants C\{uS) and C 2 (u>) 
as 



h(R) 



1 ' d*k 6 



ik-x 



(2tt) 3 J k A q 



(B.4) 



APPENDIX C: FORMULAS 



For any real positive parameter a and b, one obtains [2 

roc r^r 
/ dx x^ 1 / 2 Kq(cix) J\/ 2 {bx) = i / — sinh 
Jo V 2o 



da; £ 3/2 i^ (aa;) Ji/2^>x) = \ — 



1 



. 2 (a 2 + 6 2 ) 3 / 2 ' 

It is easy to get the following equations by differentiating Eq. (|C.l[) with respect to a, 
^ dx x^ Kl {ax)J 1/2 {bx) = ^— a{a2 + b2)1/2 , 

y o dx x5/2 ( ax ) J V2 ( bx ) = yy (q2 + b2)5/ 2 ■ 

In Eq. (|A.8p . we use the integral representation of ^0(2;) 

_ _ . . f°° , cosfzt) 



(C.l) 



(C.2) 



(C.3) 
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